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Abstract 

In this paper, we consider the interactions between a rigid body of general form 
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1 Introduction 

In this paper, we investigate the motion of a solid in an incompressible perfect fluid. The 
behavior of the fluid is described by the Euler equations. The solid we consider is a rigid 
body, whose movement consists of translation and rotation. And its motion conforms to 
the Newton's law. Assume that both the fluid and the rigid body are homogeneous. For 
simplicity of writing, the density of the fluid equals to 1. The domain occupied by the solid 
at the time is 0(t), and Q(t) = M 3 \ 0(t) is the domain occupied by the fluid. Suppose 
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0(0) = O and f2(0) = share a smooth boundary <90(or dQ). The equations modeling 
the dynamics of the system read(see also [18J ) 



— + {u-V)u + Vp = f, in 0(i) x [0,T], (1.1 

div u = 0, in x [0,T], (1.2 

u- it = (h! + u x (s-ft(t))) • n, on dd(t) x [0, T], (1.3 

lim u(x,t) = Uoo, (1.4 

m /t" = / pnda + f Tb , in [0,T], (1.5 

(Juj)' = (x-h(t)) xpnda + T rb , m[0,T], (1.6 

n(x,0)=no(x) x G $7, (1.7 

/i(0) = 0eR 3 , h'(Q) = l e M 3 , w(0) =w £ M 3 . (1.8 

In the above system, u and p are the velocity field and the pressure of the fluid respectively. 
/ is the external force field applied to the fluid. f r b and T r b denote the external force and 
the external torque for the rigid body respectively, m is the mass, and J is the inertia 
matrix moment related to the mass center of the solid. Suppose the density of the rigid 
body is p, then 



and 



m = pdx = I pdx, 
Jo(t) Jo 

[J(t)]ki = [ p[\x- h(t)\ 2 5 kl - Or - h{t)) k {x - h(t)){\ dx. 
JO(t) 

Here h{t) denotes the position of the mass center of the rigid body and 5m is the Kronecker 
symbol. And denote J(0) by J. uj(t) is the angular velocity of the rigid body, ft is the 
unit outward normal to dd(t). 

Assume that the center of O is the origin, i.e., 

y dy = 0e R 3 . 

O 

Let Q(t) be a rotation matrix associated with the angular velocity uj(t) of the rigid 
body, which is the solution of the following initial value problem: 

^ = A( W (t))Q(t) 
at (1.9) 

Q(0) = Id. 
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Here 

/ — W3 UJ2 

A(uj) = U3 — u)i 
\ -u 2 U)i 
and Id is the identity matrix. Then the domain 0{t) is given by 

0(t) = {Q(t)y + h(t):yeO}. 

For simplicity, we assume that / = 0, Uqo = 0, f r b = and T r i, = 0. 

For the case that the fluid is viscous, there have been many results over the last 
two decades. The existence of global weak solutions of the above system was proved by 
[iElElEKinilianailD]. If the rigid body is a disk in M 2 , T.Takahashi and M.Tucsnak[23] 
showed the existence and uniqueness of global strong solutions. Later, P.Cusmille and 
T.Takahashi [2] extended the result to general rigid body case in IR 2 . They also proved 
the local existence and uniqueness of strong solutions in IR 3 . 

For the case that the fluid is inviscid, [T7] dealt with this problem first. When the 
solid is of C 1 and piecewise C 2 boundary, and the fluid fills in IR 2 , a unique global classical 
solution was obtained under some assumption on the initial vorticity in [T7] . A global weak 
solution was constructed in [25] when the initial data belongs to W 1,p , p > |. Recently, 
C. Rosier and L. Rosier [TB] proved the local existence of W^ ,2 -strong solutions for d > 2, 
s > [d/2] + 2 and the solid is a ball. The key idea is to make use of the Kato-Lai theory, 
which was originated from |14| . 

At the same time, the fluid-rigid body system which occupies a bounded domain was 
also studied. In particular, [11] proved the existence and uniqueness of strong solution for 
the three dimensional case. The approach of [11] follows closely the idea in jTJ, which is 
used to study the classical Euler equations. Their method also applies to the case that 
there are several solids in the fluid. 

In this paper, we plan to extend the result of [18] to a more general setting. We will 
study the case that the solid is smooth and of a general form. As shown by the system 
itself, it is a free boundary problem. To deal with the free boundary problem, usually the 
first step is to fix the boundary. To fulfill that, [18] made a translation of the coordinate 
system. However, this special transformation can only be applied to the case that the 
solid is a ball. The solid in this paper is of general form, hence a different change of 
coordinates from |18j should be introduced. As said before, the motion of the solid is a 
rigid body movement, so a natural idea is to use a coordinate transformation consisting 
of the translation and rotation of the solid. We tried in this way. As will be shown in 
section 2, after the transformation, the new equivalent system has some term difficult to 
control. So we gave up the idea and then applied a new transformation which coincides 
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with the movement of the solid in its neighborhood and becomes identity when far away 
from it. The concrete form of the transformation will be given in section 2. In fact, this 
kind of transformation has been used by \12\ [2l [6] . 

For the new equivalent system after the transformation, we will use the Kato-Lai 
theory [H] to construct a sequence of approximate solutions and prove the limit is the 
solution required. The Kato-Lai theory is a Galerkin method in spirit. One can construct 
the Galerkin approximation by himself or herself. We use the theory here directly to avoid 
unnecessary details. By the way, the method we apply here can also be used to deal with 
the several-solids case after some minor modification. 



2 Transformed equations and main result 

As noted above, to fix the boundary, one method is to use the coordinates transformation. 
As shown in |161I17|. a direct way is the following one consisting of translation and rotation, 

x = Q(t)y + h(t), u(y, t) = Q(t) T u{Q{t)y + h(t),t), 

p(y,t)=p(Q(t)y + h(t),t), h(t)= f Q(sfh'(s)ds, 

Jo 

J = J(0), u(y,t) = Q(t) T cj, 

where Q(t) is given in section 1 and Q(t) T is the transpose of Q(t). 

After the transformation, an equivalent system is obtained as follows: 



du — 

— + [(u — h — Co x y) ■ V]u + Co x u + Vp = 0, 


rn x [0,T], 


(2.1) 


div u = 0, 


in U x [0,T], 


(2.2) 


u ■ n = (a + Hi x y) ■ n, 


on dn x [0,T], 


(2.3) 


mh" = / pnda — mCo(t) x h'(t), 
Jdii 


in [0,T], 


(2.4) 


JCd' = y x pnda + (Ju>(t)) x u>(t), 


in [0,T], 


(2.5) 


u(y,0) = n , 


y G n, 


(2.6) 


h(Q) = 0, h'(0) = l , u)(0) = wo- 




(2.7) 



The new problem is a fixed boundary problem now. However, there is a term [uj X y) ■ 
V)]n, whose coefficient become unbounded at large spatial distance. For the 2D case, the 
difficulty was overcome in |17| by assuming that uq belongs to some weighted space. The 
method there depends strongly on the fact that the vorticity satisfies a transport equation 
in 2D. However, the fact does not hold any more in 3D. To avoid this difficulty, we will 



4 



use another change of variables. The new transformation coincides with Q(t)y + hit) in 
a neighborhood of the solid and becomes identity when far away from it. In fact, this 
transformation was initially applied by Inoue and Wakimoto [12], later by Dintelmann 
etc. [6] and Cusmille and Takahashi [2] 

More precisely, for a pair of continuous vector-valued functions (l(t),Lo(t)), let 

h(t) = [ l(s)ds, t 6 [0,T], (2.8) 
Jo 

and 

V(x, t) = l(t) + u(t) x (x - h(t)), (x, t)e» 3 x [0, T], (2.9) 

which is a rigid body movement. 

Choose a smooth function £ : R 3 — > R with compact support such that £(y) = 1 in a 
neighborhood of O, and set 

= £ (Q(*) T (* - M*))) , 

where Q(t) is given by (jl.9p . 

Then introduce the functions W and A, 

1 \x-h(t)\ 2 
W(x, t) = -l(t) x(x- h(t)) + ^———ui, (2.10) 

A(x, t) = ipV + W> x W. (2.11) 
It is easy to check that A satisfies the following two lemmas (or refer to [21 [6]). 

Remark 2.1. In what follows you will see the function A will produce a transformation 
which coincides with Q(t)y + h(t) in the neighborhood of the solid and becomes identity 
when far away from it. The reason why we use A instead ofipV is that we need the function 
to be of divergence free. The function W is introduced to eliminate the divergence of i/jV. 

Lemma 2.2. (1) A(x,t) =0, if x is far away from 0(t); 

(2) A(x, t) = l(t) + x (x - h(t)) in 0(t) x [0, T]; 

(3) div A = in R 3 x [0,T]; 

(4) For all t € [0,T],A(-,i) is a C°°(R 3 ,R 3 ) function. Moreover, for every s G N, 
||A(.,t)|| W r.,cc (B a ) <C(s,T)(\l(t)\ + \ u (t)\); 

(5) For allx € R 3 , the function A(x, •) is in C°([0, T]; R 3 ) ; provided that l,u G C°[0,T]. 
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Next, consider the vector field X(y,t) which satisfies 

??M±=A(X(y,t),t), te(0,T], 
X(y,0)=y£R 3 . 



(2.12) 



Lemma 2.3. For every y G M 3 ; the initial-value problem h2.12\) admits a unique solution 
X(y,-) : [0,T] -)■ R 3 , which is a C 1 function on [0, T] . Moreover, the solution has the 
following properties: 

(1) For all t G [0,T],the mapping y i-> X(y,t) is a C°° diffeomorphism from O onto 
0(t) and from £1 onto f2(i). 

(2) Denote by Y(-,t) the inverse of X (•,£)■ Then for every x G R 3 the mapping t h-> 
Y(x,t) is C 1 -continuous and satisfies the following initial value problem, 



' dY(x,t) 

Y(x,0) = x G 



dX(Y(x,t),t) 
dy 



n -l 



A(s,t), t£(0,T], 



(2.13) 



f5) For every x,y G R 3 and /or even/ i G [0, T], the determinants of the Jacobian 
matrices Jx of X(y,t) and Jy ofY(x,t) both equal to 1, i.e., 



det (J x (y,t)) = det (J Y (x,t)) = 1. 



(2.14) 



Let 

x = X{y,t), v{y,t) = J Y (X(y,t),t)u(X(y,t),t), 

q(y,t) =p(X(y,t),t), H(t) = Q(t) T h(t), L(t) = Q(t) T l(t), 

and R(t) is the vector-valued function satisfying 

A(R(t)) = Q(t) T A(oj(t))Q(t). 



(2.15) 



(2.16) 



Denote 



G(y,t) = {g ij (y,t)) = ^ ^ X (y,t),t)^(X(y,t),t)\ , 



EdXk , s dXk , . 
^(y,t)^-( y ,t), 



dgn dgi 



2 ^ J I <% % % 



(2.17) 
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Now one can transform the original system (|l.ip - (|1.8|) into the following system, which 
is a fixed boundary problem, (see [21 [15]) 

— ^ + Mv + Nv + G ■ Vq = 0, 
at 

div v = 0, 

v(y,t)-n = (L(t) + R(t) X y) ■ H, 



infix [0, T] 



in x [0,T], 



on dtt x [0, T] 
in [0,T], 



mL'(t) = / gnder — mR(t) x L(i), 
Jan 

Ji?'(i) = / y x gnda + JR(t) x J2(t), in [0,T], 
Jan 

f(y,0) = « (y), t/efi, 
ff(0) = 0, L(0) = Iq, R(Q) = oj . 



(2.18) 

(2.19) 
(2.20) 

(2.21) 



(2.22) 

(2.23) 
(2.24) 



where 



j=i yj j,fc=i 

i=i y3 j,k=i 



i dY k dYjd 2 X k 

j ' k dt dx k dtd yj 



(2.25) 



i=i yj 

Now the main result in this paper reads 



Theorem 2.4. Suppose that s > |, no £ H s (£l) and uq ■ n = (Iq + uio x y) ■ n on dQ. 
Then there exist some To > and a unique solution (v, q, L, i?) of \2.1S\) - ^2.2J^ such that 

v € C([0,T ];iT(fi)), Vg e C([0, T ]; ff'" 1 ^)), 

and 

i.EeC'do.ro];! 3 ). 

5*1*0/1 a solution is unique up to an arbitrary function of t which may be added to q. 
Furthermore, To does not depend on s. 

Remark 2.5. Following similar proof in section 7 for the uniqueness of the solution, we 
can get the stability of the solution. The method is standard as for the classical Euler 
equations, so we omit the details. 

Remark 2.6. There are still many questions for the fluid-rigid body system. For example, 
whether the strong solution blows up or not? If there are many solids in the fluid, is there 
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any collision between these solids? Which comes first, the collision or the blow up of 
||Vu||z°= ? 

For readers' convenience, we list the outline of the following sections. Section 3 will 
be devoted to some notations or definitions for the function spaces and some lemmas to 
be used in the proof of the main result. The preliminaries contain the decomposition of 
the function spaces associated with this particular problem, the Kato-Lai theory for the 
construction of approximate solutions, and the estimates for the coordinate transforma- 
tion. In section 4, we will give the a priori estimates for the new pressure term Vg in 
(|2.18p . under the assumption that IMI#s(n)nX ^ s bounded. These estimates will be used 
for the uniform estimates for the approximate solutions, which are constructed in section 
5. Moreover, some uniform estimates for these solutions are derived at the same time. 
The convergence of these solutions is studied in section 6. And the limit is exactly the 
solution required. In the last section, section 7, we prove the uniqueness of the solution. 
Making use of the uniqueness result, the continuity of the solution with respect to time t 
in some space H s is also proved. 

3 Preliminaries 

Before stating our main result, we'd like to introduce some notations and definitions. 

Suppose 5 is a domain in M 3 . L 2 (S) is the space of L 2 -integrable functions with the 
standard inner product (•, -)l 2 (S)- By the way, we will not distinguish the scalar function 
spaces and the corresponding vector-valued function spaces. 

Suppose s is a nonnegative integer, then 

H S (S) = {u£ L 2 {S) : D a u G L 2 (S), V a G N 3 ,s.i., \a\ < a}, 
with the inner product 

(u,v) H s (s) = (D a u,D a v) L2{s) . 

\a\<s 

The homogeneous Sobolev space 

D 1 ' 2 (S) = {u G L} oc (S) : Vu G L 2 (S)}, 

with the seminorm 

MD1.2(S) = ||VU|| L 2( 5 ). 

If one identifies the two functions u±,U2 G D 1,2 (5) whenever \u\ — ^2|d 1 > 2 (5) = 0; i- e -> 
u\ and U2 differ by a constant, the quotient space D 1,2 (S), with the norm | • Ixjws), is 
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derived. In the following text, without any confusion, we do not distinguish the elements 
in D 1 ' 2 (S) and D 1 ' 2 (S) very strictly. 

Let Bn(y) denote the ball centered at y and with the radius R. := fin Br(0). Let 
P = TZTp wnere 1^1 stands for the volume of O. Hence p is the density of the solid. Let 
X = L 2 (R 3 ) be endowed with the inner product, 



X* = {u G X : div u = in R 3 , 3 i,w G M 3 , si., « = I+ wxi/inO} ) 
which is a closed subspace of X. 

Remark 3.1. For every u G X*, and suppose that u = l + uxyonO. In fact, I and oj 
are unique vectors satisfying the above relation, and 



The fact has been proved, see J3J/ or [26]. In what follows, we will denote the vectors I, oj 
associated with u G X* by l u , u u . 

Let H s = {u G X : u\q G iP(f2)} be endowed with the scalar product 

(u, v) Hs = (u, v) H s {n) + p(u, V) L 2 {0) . 

V s is the space of functions v G H s such that v\n belongs to T>(A), where A is the elliptic op- 
erator Af = Yl\ a \<s(~ l)^d 2a f with Neumann boundary conditions, and T>(A) C H 2s (£l). 
V s is endowed with the scalar product 

(u,v) Vs = (u,v) H 2s (n) + p(u,v) L 2 {0 y 

As in [18], we introduce a bilinear form onV s xX: 



Since the main idea in this paper is the Kato-Lai theory, so we'd like to give a brief 
description of the theory, which is cited from [18]. For more details, please refer to [14J. 
Let V, H, X be three real separable Banach spaces. We say that {V, H, X} is an admissible 
triplet if the following conditions hold. 




Define 
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• V C H C X, and the inclusions are dense and continuous. 

i 

• H is a Hilbert space, with inner product (•, and norm || • \\h = (•, -)jL 

• There is a continuous, nondegenerate bilinear form on V x X, denoted by (•,•), such 
that 

(v, u) = (v,u)h for all v £V,u E. H. 

Recall that the bilinear form (v, u) is continuous and nondegenerate when 

\(v, u)\ < C\\v\\v \\u\\x for some constant C > 0; (3.1) 

(v, u) = for all u£l implies v = 0; (3-2) 
{v, u) = for all v E V implies it = 0. (3-3) 

A map A : [0, T] x H — > X is said to be sequentially weakly continuous if A(t n , v n ) — 
A(t, u) in X whenever t n t and t> n — 1 t> in if. We denote by C w ([0, T]; H) the space of 
sequentially weakly continuous functions from [0, T] to H, and by C^([0,T]; X) the space 
of the functions u G VF 1 ' oo (0,T; X) such that ^ G C w ([0,T],X). 

We are concerned with the Cauchy problem 

dv 

— + A(t,v) = 0, with v(0) = vq. (3.4) 
at 

The Kato-Lai existence result for abstract evolution equations is as follows. 

Theorem 3.2. Let {V, H, X} be an admissible triplet. Let A be a sequentially weakly 
continuous map from [0, T] x H into X such that 

(v,A{t,v))>-p(\\v\\ 2 H ) for te[0,T], v G V, 

where j3{r) > is a continuous nondecreasing function of r > 0. For any vq G H , consider 
the ODE j'(t) = /3(7(t)), 7(0) = ||vo||/fv Suppose the maximal time of existence 0/7 is 
Tq, then for \3. 4\ ), there exists a solution v of fi3.4\ ) in the class 

veC w ([0,T o ];H)nCi([0,T o };X). 

Moreover, one has 

\\v(t)\\ 2 H<l(t), te[0,T o ], 

In fact, it was proved in [18J that the triplet {X , H S ,V S } is admissible for any smooth 
boundary. 

The following lemma gives a decomposition of L 2 (IR 3 ). In particular, the second part 
is a replacement of Proposition 3.1 of |18j . 
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Lemma 3.3. Let G\ = {uG L 2 (M 3 ) : u = Vq x , qi G L^QR 3 )}, and 

G\ = {u G L 2 (R 3 ) : div u = in R 3 , u = Vq 2 in Q, q 2 G Lj oc (Cl), u = <f> in O, 



(fidy = — I q2n da, and / (fix ydy = — I q^n x y da 
o J an Jo J an 

Then (1) X*,G\ and G\ are mutually orthogonal in the sense of the standard L 2 -inner 
product and 

L 2 (R 3 ) = X* ® Gl (B G 2 2 . 

It means that for every u G L 2 ' 



J ui + Vgi + Vg 2 , y G - 2 2 

= < G^eGjeGj. (3.5) 

[ ui + Vgi +<p, y GO 
Suppose u\ = l Ul + w ui x y in O, then 

\lui \ + l^uj < CII^iIIl 2 ^) < C||n|| i 2( K 3) < C||u||^. (3.6) 



(2) Denote the projector which maps L 2 {R 3 ) to X* by P. In fact, P maps H s into H s 
continuously for any s > 0. 

Proof. The orthogonal decomposition in (1) has been proved in [3J. And the estimate 
(|3.6p is partly derived in [3] and partly due to the fact that || • ||£,2(r3) and || • ||^ are 
equivalent norms. Now we verify that (2) holds. For every u G H s (s > 0), suppose that 
u = u\ + Vgi + V02 = u\ + VP over Q, then it suffices to prove that 

||VP|| H . ( n) < C\\u\\ H ., 

with some constant C independent of u. 

In fact, from the formula (|3.5p . P satisfies the following equations: 

AP = div u, in 

dP , , ( 3 - 7 ) 

— = u ■ n - (l Ul + uj Ui xy) -n, on dVL. 
on 



Let 



ip = V x 



• (^i x y - w «ily| 2 ) 



2 

where £ is a cut-off function defined in section 2. Clearly, div ip = in Q and (p ■ n = 
(l Ul + w ui x y) • n on <9f2. Therefore, (13. 7|) can be rewritten, 

AP = div (u — <p), in 

dP ( 3 - 8 ) 
— : = (u — (p) ■ n, on dfl. 
on 
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The solution to the system (|3.8|) is closely related to the Helmholtz-Weyl decomposi- 
tion. As proved in [TJ, 



l|VP|| H . (n ) < C\\u - tp\\ H »(n) 

< C (Ikll^(Q) + IMI#«(n)) 
<C (\\u\\ HS(n) + |Z U1 | + \u Ul \) 

<C\\u\\ H ., 

which completes the proof of Lemma 3.3. □ 

The following lemma is to give the bounds of the terms which appear in the system 
(|2.18p - (|2.24p . Before that we'd like to give another description of relationship between 
(7,a;) and (L,R), which is different from (|2.15p - (|2.16p . Suppose (L,R) is given, we want 
to determine (I, to) and then define the other coefficients in (|2,18p - (j2.24p . 

Since 

dQ(t) 
dt 

multiplying by Q T (t), then 



A(u(t))Q(t), 



dQT(t) Q(t) = A(R(t)). 



dt 

It gives that 

{¥ = Q(t)A(R(t)), 
\ Q(t) = Id. 

Now we see that if R(t) is given, Q(t) can be determined. Then l(t),u>(t) are determined, 

l(t) = Q(t)L(t), A{uo{t)) = Q(t)A(R(t))Q T (t). (3.9) 

Then A,X,Y,gij,g i:i ,T,Mv,Nv can be determined as in section 2. 

Moreover, if L(t), R(t) G C[0,T], then l(t), u)(t) £ C[0,T], with the estimate 

\l(t)\ + \co(t)\ < C(T) (\\L\\ Lao( p tT) + \\R\\l°°(o,T)) • (3.10) 

Suppose that (Z 1 ^),^ 1 ^)) and (I 2 (i) , u? (t)) are determined by (L 1 (t) , R 1 (t)) and (L 2 (t) , R 2 (t)) 
respectively in the above way, then 

\\h — h\\L°°(0,T) + 11^1 ~~ ^2||l°°(0,T) 

< C(T)(1 + \\Rl(t) ||l°°(o,T) + 11-^2 (*)||z<»(Q,T)) ■ (ll-^i - L 2 \\l°°(0,T) + \\Ri - R2\\l°°(0,T))- 

(3.11) 
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Lemma 3.4. Assume that v is a function in L°°(0,T; X*) and s is a nonnegative integer. 
Suppose there exists M* > 0, such that \\v\\ Loa r QT . x \ < M*. Let 

L(t) = l v (t), R(t) = U v (ty 

Suppose l(t),u(t) is given by L(t),R(t) as in \3. 9\) and A, X,Y, g^ , g l i ,T are defined as in 
section 2. Then for every t £ [0, T], the following estimates hold: 

\\Jx(;t)\\w:°°(a*) < C(s,T,M,), \\J Y (X(;t),t)\\ We , aom < C(s,T,M,), (3.12) 

\\A(X(-,t),t)\\w^m ^ C(s,T,M*), ||^(-,t)|| W r.,c 0(R 3 ) < C(s,T,M*), (3.13) 

\\g ij (;t)\\w^m < C{s,T,M*), \\G-\-,t)\\ ws ,oo m < C(s,T,M*), (3.14) 

||r(-,t) \\ W s^ m < C(s,T,M*), (3.15) 

where G~ 1 is the inverse of G. 

Proof. For every j = 1,2 or 3, let z(y, t) = 

dz(y,t) dA.^ r . . . . . 
-^l = -(X(y,t),t). Z (y,t), 

z(y,0) = ej, 

where is the j-th vector of the basis of M 3 . Then 

f f dA 

z{y,t) = ej + J — (X(y,t),t)-z(y,s)ds. (3.16) 



It follows from Gronwall's inequality that \z(y, t)\ < C(T, M*). 

Since det(Jx) = 1, then Jy = (J^), where J x is the cofactor of Jx- Hence 

\J Y (X(;t),t)\ <C(T,M»). 

Furthermore, 

\G(-,t)\ < C(T,M*), \G-\;t)\ < C(T,M*). 
Denote D@ = J^. From (I3.16p . one can deduce that 

Following the preceding process, one can get the estimates (I3.12p - (l3.15p . □ 
Next lemma is about the Lipschitz continuity of the coefficients with respect to v. 
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Lemma 3.5. Assume that the assumptions of Lemma 3.4 hold for v % , i = 1, 2. Let L l (t) = 
l v iM and R l (t) = 0J v iu\. Define Q l (t) , l l (t) , uJ 1 (t) and other terms correspondingly. Let 
L(t) = L l (t) - L 2 (t), R(t) = R^t) - R 2 (t), X = X 1 - X 2 , A(y,t) = A{X 1 (y,t),t) - 
A(X 2 (y,t),t), G = tfi) = (gM-g**), 9ij = g^-gl, G' 1 = (G 1 )' 1 - (G 2 )-\ and 

T L,k = T m,k - r m 2 fc' Then f° r eVer V f G [0. T }> 

\\X(-,t)\\ W s^ m < C(s,T,AQ (}\L\\ LOB(0)T) + \\R\\l°°(o,t)) , (3.17) 

||A(., t)\\ ws ^ m < C(s,T,M*) (\\L\\ L ^ T) + p|| LOO(0 ,T)) , (3-18) 

||^(-,i)|| Ws ,°o (K3) < C(s,T,M*) [\\L\\ La o i0jT) + ||Ji||x=c (0)r) ) , (3.19) 

\\g l i(-,t)\\ W s,o« m < C(s,T,M*) (\\L\\ Laa{0>T) + \\R\\ L ~ {0 ,T)) , (3-20) 

WG-^QWw^ps) < C(s,T,M*) (||L|| L ao (0>T ) + \\R\\ L o° (0 , T) ) , (3.21) 

|[r(-,*)|]wr-.~CR 3 ) ^ C ( S ^ M *) {\\L\\l°°(o,t) + \\R\\l°°{o,t)) • (3.22) 

Proof. Since 

9 -X^ll=K(X\y,t),t)-k{X 2 {y,t),tl 
X(y,0) = 0, 

Simple calculation and the estimates in Lemma 3.4 induce that 

\A(X\y, t),t) - A(X 2 (y, t),t)\ < C(T, M^X^y, t) - X 2 (y, t)\ 

+C(T,M*)(\L(t)\ + \R(t)\). 

Therefore, 

II^I|C([0,T];L°°(K 3 )) < C(T,M*) (||L|| LO o( 0jT ) + ||-R||i°°(0,T)) • 

\\My,t)\\Lo°(W) <C(T,M,)\\X 1 (y,t)-X 2 (y,t)\\ L ^ m 
<C(T,M*)(\\L\\ Lao{0tT) + \\R\\ La . {0tT) ). 
Other estimates can be derived similarly. □ 



4 A Priori H s - estimates of 

In the following text, s > |. Given a function v S C w (0, T; H s n X*), which satisfies that 
IMlL°°(o,T;.ff s ) — Mq. Suppose v is a solution to (|2.18p - (|2.25p . taking the divergence of Eq. 
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(|2.18p . then q satisfies the following system, 
-,dq_ 



div I ^2g lj p~ I = -dw(Mv + Nv), in O, 



V g lj p-m+ (— [ qnda) ■ 



n + ( J I y x gnd<7 ) x y • n 



= —(Mv + iVv) • n + (w„ x / tI ) • n — [ J 1 (Jw„ x u v )] x y • n, on 50. 

Here Q, Z, co, g ij , and iVu are given as in $$3$, ([2Tl1^ - (j2TTTl) and (^25]) . replacing L, i? 
by Z„,u;„. 

In this section, we will give the IP-estimates of Vq, for every fixed time t E [0, T]. For 
simplicity of writing, we omit t. Here is the main result of this section. 

Proposition 4.1. 

l|V?|| H .(n) <C(T,M )(l+\\v\\ Hs ), 
with some constant C depending on T and Mq. 

Remark 4.2. The system q satisfies is almost a classical elliptic problem with Neumann 
boundary condition. During the proof of Proposition J^.l, the main idea is to use the 
Lax-Milgram theorem to get the existence of weak solution, and the standard high-order 
regularity estimates for exterior elliptic problems. 

Proof. For every fixed t 6 [0,T], the matrix G = (g^) = JyJy, so G is positive definite. 
Denote Aj(y, t) > 0, (i = 1, 2, 3) the eigenvalues of the matrix (g lJ ). Since det(g l:] ) = 1, thus 

3 3 3.. 

lli=i = 1 and E = E ff" > °- Let 70 = sup \g u \, then we have 370 > \i > v^-vj for 

i=i i=l y eR 3 (1 ' 

every i = 1, 2, 3. By virtue of Lemma 3.4, there exist constants C\(T, Mq) and Ci(T, Mq), 

C 1 (T,M ) < \\ i \<C 2 (T,M ), i = 1,2,3. 

Next, we shall use the Lax-Milgram theorem to prove the existence of the solution of 
(|4.1|) . and then give the iP-estimate of this solution. 

Set a bilinear form B and a linear functional F on D 1,2 (Q) as follows, for every n,q £ 
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L 2 {n) 



i=i \j=i 
+ ( •/ 



£ndcr I 
an ) 



y x rjnda I , 



y x qnda 
an / \Jan 

F(rj) = - / (Mo + iVu) • V7?dy + / (ca„ x • ^ndcr 
./n Jan 

[J (Ju^ x x y • rjnda. 

an 

Note that 

"<«. « ) = e(e + ^ (I 2 + J -w ' »<«> 

i_1 V- 1 / (ci) 

>C 1 (T,M )\\Vq\\ 2 L 2 m + -( [ qnda) + Jw(q) ■ w(q), 



(4.2) 



™ \J 9 n 



where w(q) = J 1 J 9n y x qnda. 

Since J is a positive definite matrix, then there exists some constant a > such that 

a^ 1 \w(q)\ 2 < Jw(q) ■ w(q) < a\w(q)\ 2 . 

Combining the above inequality and (|4.2p . one gets that B is coercive. 
On the other hand, 



L 2 (n) 



< II^Hl^ (R3) || Vg|| L 2 (n) || VryU^^). 



Then along the line of the proposition 3.3.1 in |25| . one can easily verify that the bilinear 
form B is bounded. 

Now we turn to the functional F. 



(Mv + Nv) ■ Vr]dx 



< \\Mv + Nv\\ L 2 {n) • ||V77|| L 2( n) 

< C (\\A\\ W i ia0 (n) + ||^|U°°(n) + ll^-x||L°°(n) + IMIl°°(c) 

+ ll r llL-(n)) • \\v\\m(ci) ■ l|Vr?|| L 2 (n) 
<C(T,M )\\V V \\ L 2 {n) , 



(4.3) 
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Choosing some r\ £ D 1,2 (Q) such that f n rjdy = 0, for some large r. 



{uj v x l v — [J (Joj v x uj v )\ x y} • r/ndo" 



< C(|w v ||Z„| + \u v \ 2 )\\ri\\ L 2 {dn) 

< C(n)(\u) v \\l v \ + \u} v \ 2 )\\r]\\ H i {nr) 
<C(T,M ,n)\\Vv\\ L 2(a), 



(4.4) 



where the last inequality is given by the Poincare's inequality. 

From the above estimates, it follows that F is bounded. According to Lax-Milgram 
Theorem, there exists a unique q € D 1,2 (£l) such that 



with 



B(q, V ) = F(ri), Vijel) 1 - 2 ^) 
\\Vq\\ LH n) <C(T, Mo). 



(4.5) 



Let 



Then 



1 f 

L\ = — / qnda, w = J 



1 1 y x qnda. 

an 



\Li\ < C||Vg||£2(n) < O(T,M ), \w\ < C\\Vq\\ L , m < C(T,M ). (4.6) 

Now we go to the .£P-estimate of Vg. Similar to [18J, the method is a standard 
regularity estimate for an exterior problem of elliptic equations. Consider the Neumann 
system which is equivalent to (14. ip . 



div 



dq 



-div(Mv + iVv), 



in £1, 



3 



<9g 



rii = — {Mv + iVf ) • n — (u v X l v ) ■ n 



(4.7) 



+ [J" 1 (Jo;„ x u;„)] x y ■ n — L\ ■ n — (w x y) ■ 



n, on 



90. 



To estimate || Vg||#s(Q), the key is to estimate the terms ||div(Mu + Nv)\\ H s-i^ and 
- (Mv + Nv) - nil i 
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< 



^dyAdt +Vl ) dy 3 

d 2 Yj dX k dvi 
i fj? =1 dtdx k dyi dyj 



+ 



•s—^ dvj dvi 



<C (ll A lliy s >°°(R3) ll^xllvy— i,^(R3) ||V«|| ff »-i( n ) + || Vu|| ff »-i(n) 11^11^0(0)1 



and 



E 

i=l 



dyi 



E^ 



dE fe , dY k 8 2 X k \ 



j ' k dt + dx k dtdy 



J i,fc=i 



<c[( 



W/s.°°(R 3 ) + II llw s >°°(R 3 ) ll"^llw s >°°(IR 3 )) ' ll w llfl-*(n) 



+ liri 



Mlff s (fi)IMIff s (n)] • 



Hence, 



Denote 



\\6to(Mv + Nv)\\ H .- Ha) <C(T,Mo)\\v\\H.- 



\j,k=l 1 



dY k , SE fe 3 2 X fc l 



+ 



Then 



\ 

-(Mv + iVu ) • n = ( + u J • Vw • n + / • n. 
Note that there is no derivative of t> in the term /, so it is easy to handle. 



< C [(Il-T|| iy s -°°(R 3 ) l|-^-|lvy s '°°(IR 3 ) + ll^y lliy s ' 00 (R 3 )ll A -llw s '° (M 3 )) 

+■ ||r 

<<7(r,M Q )|Hk. 



I^IIh^o) + l|r||w s >°°(R 3 ) ll u llif=(n) ll u llfl-»(n) 



To estimate || (§E + v) ■ Vv ■ n\\ „ 1 



we proceed as in [HQS]. Note (^jr+v)' 
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on <9fL then 



9Y \ r, - 

— — h f • Vf • n 
at J 

dY \ ^fdY \ _ (BY \ „(dY\ _ 

it +v )- v {-w +v )- n -{-w +v )- v {-w)- n 

E (jPl + " (Z + • V ■ ft. 



Hence, 



dY 
~dt 



+ v } ■ Vv ■ ft 



< C 



Combining (|4.9p and (|4.10p . one has 



2 



+ l) (1 + \\v\\ H . (a) ) 



(Mv + Nv) ■ n\\ H s-^ m) < C(T, M )(l + \\v\\ Hs ). 



The other terms can be estimated as follows: 



[J-^Juxu)] Xyfi\\ HS _i m <C(n,R)\u;\ 2 <C(T,M ), 



\\L,-n\\ H ^ m) <C\L l \<C{T,M ), 



\\{u v x l v ) ■ A\ H s-\ m < C(n,R)\u\\l\ < C(T,M ), 



IK™ x y) -n\\ H s^ m < C(n,R)\w\ < C(T,M„). 
Choose some r > such that supp(£) C Br, and a cut-off function £i, 

6(y) = < 



1, if |y| < 2r, 
0, if |y| > 3r. 



(4.10) 

(4.11) 

(4.12) 

(4.13) 
(4.14) 
(4.15) 
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Hence, p\ = t^\q solves the following equation 



,dpi 



^ \dyj dyi dyi dyj dy, L dyj 



in B Ar \ O, 



E 



E 



•(¥« + JVv) ■ n — u> v x l v ■ ft 



+ [J 1 (Jlj v x w v )] x y ■ n — Li ■ n — w x y ■ ft, on 90, 



dpi 



rii = 0, 



on 3Ba t . 



(4.16) 



By virtue of the regularity theory for elliptic equations [22], for any (3 > 1, 



IMIi^+^B^O) ^ hi (||G|| WSj cx=( R 3), (370)^ 



E 



^-l(B4r\0) 



||eidiv(Mt; + iVt;)||^- 1(B4AO) + 



y> ij ( dj\ dq d£i dq 
\dyjd yi d yi d yj 



HP-^{B ir \0) 



(4.17) 



d 2 i 



dyidyj 



+ || - (Mv + Nv) ■ n\ 



HP-^(B ir \0) 

+ \\J- x (Juj x uj) ■ n\\ Hf3 _i {aa) + \\{Lx +w x y) • n||^_i (an) + \\pi\\mB 4r \o) y 
where hi(-, ■) can be chosen an increasing function with respect to both variables. In fact, 



E 



dg ij dft 



i-j-- 



\ dyi dyj 



+ 



H0-l(B4r\O) 



i;j 1 d£i dq | dCi dq 
^ ' dyj dyi dyi dyj 



HP-l(B4r\0) 



+ 



d 2 i 



Y * 

£+y dyid yj 



(4.18) 



W-i{B Ar \0) 

< C (\\G\\ W s,oo(^i) ■ \\q\\Hf ) - 1 (B 3r \B 2 r) + l|G ! |lwP- 1 . 00 CR 3 ) ' ll^||^-l(B 3 ».\B2r)) • 

Combining the above estimates, one gets 

1 (B 2r \0) 

<C(T,M ,r) (||g||^-i (B3AB2r) + \\q\\ L HB 3r \0) + l) • 



(4.19) 
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Choose some particular q such that 

q(y)dy = 0. 



'B 4r \0 

It is reasonable, since q is still a solution to (|4.7p if it is added by a constant. By Poincare's 
inequality, 

IMIf^+i(B 2r ) ^ C(r,M ,r) (||g||^/3 ( B 3r \ij 2r ) + ||Vg|| i2(B4 x C) + l) 

• " v 7 (4-20) 

<C7(r ) Mo J r)(|| g [|^(fl Sr \B 3P ) + l)- 

It implies that high-order regularity of q can be controlled by the lower-order regularity. 
Therefore, using this method by choosing appropriate r, we can get that for every R > 0, 
such that O C Br , 

\\Vq\\ H s(Si R )<C{T,M ,R)(l + \\v\\ Hs ). (4.21) 
Fix some R big enough. Choose some smooth cut-off function £2; such that 



0, if \y\ < lR, 

1, if \y\ > R. 



6 GO = 

Since q %] = 5^ outside Br , hence 

J 2 

div(G • Vq) = Aq. 

Let p 2 = then 

Ap 2 = 6("div(Mw + Nv)) + 2V£ 2 • Vq + A&q := /. (4.22) 

Therefore, 

l|Vp 2 ||^(R3) < C (l|/l|^-l(K 3 ) + I|VP2||l 2 (K 3 )) • ( 4 - 23 ) 

/ is estimated as follows, 

< C (\\div(Mv + Nv)\\ H .- Hn) + || ^« II ja-— x ( « < Jvl + <!„[<*)) ( 4 - 24 ) 

<C(T,M )(l + \\v\\ Hs ). 
Hence 

||Vg||^3 (M 3 XjBR) < C(T,M )(1 + \\v\\ Hs ). (4.25) 
(|4T211 and (|4T25l) give that 

||Vg|| H . ( n) < C(T,M )(l + \\v\\ Hs ). (4.26) 

It completes the proof of Proposition 4.1. □ 
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5 Construction of approximate solutions 



In this section, we will construct a sequence of approximate solutions. Similar to [18j . the 
main idea is the Kato-Lai theory. Suppose for every n > 0, there is a function v n {t) € 
C w ([0,T];H s ). Denote L n (t) = l ¥v n (t) , R n (t) = u Pv n (t) . Then L n (t),R n (t) € C[0,T\. 
Solving the following initial value problem 



dQ n (t) 



dt 

Q n (0) = Id. 



Q n (t)A(R n (t)), 



(5.1) 



One can get a solution Q n (t). 
Define 



r(t) = Q n (t)L n (t), A(u; n (t)) = Q n (t)A(R n (t))[Q n (t)] T 
r = i{{Q n ) T it){x-h n {t))), h n = f Q n {s)l n \s)ds 



and 



IT" 



n in 



V H = l n (t)+LO n (t) X ( X -h n (t)), 

l n {t) x (x - h n {t)) + l - -^-u n (t), 



where £ is a cut-off function given in section 2. 

Let A n = i) n V n + Vip n W n . Hence one can define X n (-,t), A n (X n (-,t),t), gV' n , g%, 
rj.£, M n ,N n given in OTTl) and (l235|) . Suppose that q n is the solution to the following 
system, 



ZpL d yj m Van J 



div(M n ¥v n + N n ¥v n 

J 1 I V x q n nda 



in $7 



mi 



x y ■ n 



-(M n Pv n + N n Fv n ) -n-(R n xL n )-n + [,J- l (JR n x R n )] x y ■ n, on 

(5.2) 
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Now define an operator A n (t,v) as in 



A n (t,v) = l Q 



at 



+ Fv n ~ 1 ) • Vu- Q 



In 



dY 



n—l 



+ Fv™' 1 ) • VFv 



dt 



+ . 



me > 

j,k=i i 



dt dxk dtdyj 



r 1 + E r ;: 



r\wv)]-\Fvr k -' 



n-l 



+ 



"» Jen 



j,k=i 

+ IP 1 (i?™- 1 x L n_1 - J- l (JR n - 1 x i?"- 1 ) x y) 



lo 



nda — [ J 1 / y x q n 1 nda I x y 



(5.3) 



where the operator Q = I 
Let 



Mv) 



m(y), y g fi, 



Consider the following Cauchy problem, 



u? + A n (t,i; n ) = 0, 



(5.4) 



where v o € H s n X* . 

In particular, let v°(y,t) = Vo(y). We shall prove that for each n € N, there exists a 
solution v n £ C w (0,T n ; H s ) with some uniform lifespan T n . 
For simplicity, denote 



(F n -^) 



+ Pt."- 1 ] • Vvi, (G n ~ l ■ Vq n - l )i = V ^' n - 1 4 • 



dt 



En 



dt dxk dtdyj 



,.n-X\ 



+ E rg- 1 ^- 1 ^^ 





f 1 


lo 










" 1 




m 



yn— 1 .. rn—1 t— 1 / rryn—X nn—X\ 



„n-l. 



nda — I (J y x q n nda) x y 
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and denote M = 2\\vq\\h s - Suppose that there exists some To > such that for all k < n, 

\\v k \\L°°(p,T ;H a ) < M o- 

For the estimate of (v, A n (t, v ))h s , 

\(v,A n (t,v)) Hs \ < \(l n F n - 1 v,v) Hs \ + \(Q(l n F n - 1 Fv),v) Hs \ + ftFET^v, v) Hs \ 
+ |(P(l fi G"- 1 V 9 "- 1 ),^) ff | + \(FK n - 1 ,v) Hs \ 
:= J1 + J2 + J3 + J4 + h 

Then we estimate them term by term. Starting from the easiest one, 
J 5 < \\¥K n - 1 \\ Hs \\v\\ Hs 



< c 



q n ~ l nda 



an 



+ 



< CdL"- 1 ] • \R n -\ 
<C(T ,M )\\v\\ Hs . 

By Lemma 3.3 and Lemma 3.4, 



y x q n 1 nda 
an 

j^n—l\2 1 ||Y7„n— 1 



||V<f 



+ il"- 1 !!^"- 1 ! + li?"- 1 ! 2 ^ ||«|| Hi (5.5) 

L2(n)) • M\h 3 



J4<||P(ln(G"- i Vg"- i ))||^-||t;||^ 



< C||G n-1 Vg n-1 | 



H*(n) ■ \\ v \\h s 



(5.6) 



<C(T ,M )\\v\\ Hs . 
J3 is also easy to estimate since there is no derivative of v or w n_1 , 

h <C(T ,M )\\v\\ Hs . (5.7) 
Now the most difficult terms J\ and J2 are left, since there is derivative of w or Pu. 

3 /> / o-i^n— 1 



\ct\<s ct\<a i=l 

When a\ = (0,0,0), since 



dY 



+ Fv n ~ 1 ) ■ Vd a ~ ai Vid a Vidy 



dt 



(5.8) 



div 



dt 



n-l 



in CI, and 



dY 



n-l 



dt 



+ ¥v n ~ L • n = on dft, 



then 



ay 



n-l 



+ Pv n ~ 1 ) ■ Vd a Vid a Vidy = 0. 



dt 
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Therefore, we assume that \a\\ > 1. Let Q2 = a — a\ 
■+Pv n - 1 ) -Vd a2 Vi 

• Vd a2 v 



< 



301 



at 



L 2 (n) 



L 2 (Q) 

+ ||0 a W- 1 -V0 aa w i || La(n) 



< C{T Q ,M )\\v\\ Hs + ||^ n_1 ||Loc (n) ||Vx;||^-i (Q) + IIVwIIlco^IIP^Hh.o 
<C(T ,M )\\v\\ Hs . 

Hence, 



Ji <C(T ,M )\\vf Hs 



For the term J2, 



J 2 < ||Q[r n F n - 1 Pu]||H.(n)ll«llH.(n) + C^-^^M^y 



Herein, 



ay-n-1 



at 



<C(T ,M )\\v\\ 2 Hs . 
To estimate || Q[1^^- 1 P^] || ijs(H), consider the following system 



MI#JMIl 2 (0) 



A(f> = div I — + Pi;"- 1 J • VPu, in fi, 

a^A /ay n_1 \ 

^- = + ¥v n - 1 J • VPv ■ n. on 

on \ at J 



In fact, Q^n-P^Pu] = V0. 
Note that in the domain fi, 



div 



dY 



at 



+ Pf™ -1 ) • VPv 



i.i=l 



-^—+Pv) di(Pv)j, 



and on the boundary a$7, 

+ P^" 1 ) • VPt; • n = + Pi;"" 1 ) • V (P« - ^ - x y) 

+ (^T 1 + Pun_1 ) • V + wp, x y) • n 
Hence, as estimating Vg in section 4, one can get that 

\\V<t>\\ H s(n)<C(T ,M )\\v\\ Hs , 
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consequently, 

J 2 <C(T ,M Q )\\vf Hs . 

Therefore, we have 

\(v,A n (t,v)) Hs \<C(T ,M )(l + \\v\\ 2 Hs ). (5.10) 
Now, fix some big To. Let's consider the corresponding ordinary differential equation, 

7 '(i) = C(r ,Afo)(l + 7(t)), 7(0) = \\vof Hs , 

with C(Tq,Mq) the same constant as in (|5.10p . Assume that T < Tq is a time such that 
for every t € [0, T], 

7(t)<4\\v \\ 2 Hs = (M ) 2 . 

Then by the Kato-Lai theory, the solution v n to (|5.4|) can be derived at least on [0,T], 
i.e., v n € C w {[0,T];H s )nCU[0,T];X), and 

\\v n (t)\\ 2 Hs <l(t)<(M ) 2 , te[0,T}. 

For n = 1, we choose v°(y,t) = vo(y). Following the above process, one can construct 
a solution v 1 to the system (|5.4p . By iterating the same steps, a sequence of approximate 
solutions {v n } can be constructed. 

6 The convergence of approximate solutions 

In this section, we show that {v n } converges to a solution of the system ()2.18j) - ()2.24p . 
According to the estimates in section 5, 

IKIIl~(o,iw s ) < M , (6.1) 

H^1l~(o,t ; x) <Mi, (6.2) 
l|Vg n ||Loo( 0jT; ^-i(^)) < M 2 . (6.3) 
Since P is a bounded operator on H s , X, and it commutes with dt, then 

l|^ n ||L-(o,T ; if s ) < M 3 , (6.4) 
||^P^|| ioo(0)T;1) < M 4 . (6.5) 
Hence from the [21], there exists some function v G C w ([0,T]; H s ) such that, 

v n in C W (Q,T;H S ), (6.6) 
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Fv n ^Fv in C w (0,T;H s ). 
By the Aubin-Lions lemma, for every ro large enough, 

v n ^v in C{[Q,T];H s - l {Sl ro )rM?{B ro )), 

Fv n -^-Fv in CdO^li^^^nL 2 ^)). 
Moreover, (|6.9p implies that 

L n (i) -> L(t) = Z Pu(t) in C[0,T], 

R n (t)^R(t)=u ¥v{t) in C[0,T]. 
While (|6.3|) . (|6.8|) and (|6.9|) tell that there exists some function q such that 

Vq n ^Vq in C w ([0,T];H s ~ l (p.)), 



q n nda 



Ml 



qnda in C[0,T], 



/ y x q n nda ->• / y x gndcr in C[0,T]. 
Van Van 



(6.7) 

(6.8) 
(6.9) 

(6.10) 
(6.11) 

(6.12) 
(6.13) 

(6.14) 



In fact, g is a solution to the system (|4.1|) . It can be seen by taking the limit of (|5.2I) , 
From all the convergence results (|6.ip - (l6.14p . it follows that 



v t + A(t,v) = 0, 
v(0) = v , 



(6.15) 



where 



dY 



A(t,v) = 1 Q [—+Fv)-Vv-Q 




} (p u ) i + ^ r 



J" 



l'o (RxL + J-\jR x R)xy) 
y x qnda I x y 



Next, we shall prove that f is a solution of the systems (|2.18p - (|2.24p . The proof starts 
with the observation that v(t) = Fv(t), for all t £ [0, T]. In fact, applying Q to each term 
in (|6.15p and taking the inner product with Qv(t) in X yields 



--\\Qv 



x 



+ (Qv(t),QA(t,v))x=0. 



(6.16) 
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Note that div + Fv) = in U and + Fv) ■ n = on <90, then 
(Qv(t),QA(t,v))x = (ov(t), In (^+^) -VOu® 



(6.17) 



Since vq = Fvq, it tells that Qvq = 0. Hence, for every t G [0, T], Q"u(i) = 0. Therefore, 
(|6.15p can be written as 



dv -> 

— + P[l n (Mw + Nv + G ■ Vq)) + 



I'd ( — R x L - J- l UR x R)xy 
m 



(6.18) 



Jan 



qnda — (J / y x gndo") x y 



no 



Taking the inner product in X with a test function <f> G JT*, one has 

/ (i/ + Mti + Nv + G • Vq) ■ <pdy 
Jn 

+ mL' ■ l<p — qnda ■ 1$ + m(R x L) ■ 1$ 
Jan 



JR' J[J~ l (JRx R)]- U4l - J J -1 / ; 



y x qnda 



(6.19) 



0. 



For every function 4> G Cq° (K 3 ), with supp(<p) C $7, and div = in K 3 , (|6.19p yields 

/ (V + Mi) + Nv + G • Vg) • <pdy = 0, 
Jn 

After the theory of Helmholtz-Weyl decomposition, there exists a function p such that 
Vpe^^T;^- 1 ^)) and 

w' + Mw + iVw + G- Vg + Vp = in U X [0,T]. (6.20) 

From the identification of q and (|6.20p . one knows that for every t G [0, T], 

Ap = 0, in f2, 

^ = 0, on dVL. 
on 
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The above system has only constant solutions, thus 

v' + Mv + Nv + G- Vq = in Q x [0,T]. (6.21) 
Now taking some test function 4>{x) E X such that cp(y) = in O, then 

mL' ■ — [ I qnda I • + (mR x L) ■ = 0. 

\Jdn J 

Since is arbitrary, then 

mL' = / qnda - Rx L. (6.22) 
Jan 

Similarly, taking some test function (p(y) 6 X such that cp(y) = oj^ x y in 0, then 
XR' • - J{J~ 1 {JR x R)) ■ - J ( J' 1 yx qnda) ■ ui^ = 0. 

Thus 

JR' = I yx qnda +JRxR. (6.23) 
Therefore, (v,q,L(t),R(t)) is a solution to the system ([2TT8l) - (pT2l ). 

7 Uniqueness and continuity with respect to time 

In this section, we will prove that the solution of the system (|2.18p - (|2.24p is unique and 
then get the continuity in H s with respect to time. 

Assume that there exist two solutions v , v 2 € C w ([0, T]; H s ) n C^([0,T]; X) to the 
system (l2TT8j) - ([2T2i]l . then 

■Uj 1 + M 1 ?; 1 + A 1 ?; 1 + G 1 V^ 1 = 0, inOx[0,r], (7.1) 

v 2 t + MV + AV + G 2 Vg 2 = 0, in Q x [0, T]. (7.2) 

Let H 1 = (G 1 )- 1 and H 2 = (G 2 )" 1 . Multiplying ([71]) and £l2D by tf 1 and H 2 
respectively, and denote K = max{|jv 1 \\l°°(o,T;H s )-> ll^ 2 llL oo (o,T;_ff s )}- 

Subtracting the two equations and taking inner product in L 2 (VL) with v 1 — v 2 , then 
one gets 

= {H x vl - H 2 vlv l - v 2 ) L 2 (n) + (Vg 1 - Vq 2 , v 1 - v 2 ) L 2 {Q) 
+ {H l {M l v 1 + A 1 !; 1 ) - H 2 (M 2 v 2 + AV), v 1 - v 2 ) L 2 (n) 
■-h + h + h, 
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(7.3) 



Denote l v i , oo v i , l v 2 , u v 2 by L 1 , R 1 , L 2 , R 2 respectively. 

h = I V^ 1 - q 2 ) ■ (v l - v 2 )dy 
Jn 

= [ (q 1 - Q 2 )^ 1 - v 2 ) ■ nda 
Jan 

= [ (q 1 - q 2 )^ 1 - L 2 ) ■ nda + [ (q 1 - q 2 )(R} - R 2 ) x y ■ nda 
Jdn JdQ 

= m{L x - L 2 )' ■ (L 1 - L 2 ) + m{R l x L 1 - R 2 x L 2 ) ■ (L 1 - L 2 ) 
+ J(R l - R 2 )' ■ (R 1 - R 2 ) - (JR 1 x R 1 - JR 2 x R 2 ) ■ (R 1 - R 2 ) 

= ^-^ + i£ [<*#-#))■<*-*»)] 

+ miR 1 x L 1 - R 2 x L 2 ) • (L 1 - L 2 ) - (JR 1 x R 1 - JR 2 x R 2 ) ■ (R l - R 2 ). 
The term I\ can be estimated as follows, 
h = {H l v]-H 2 v 2 y-v 2 ) L , m 

= (H l {v l -v 2 ) u v l -v 2 ) L2m + ((H l -H 2 )v 2 y-v 2 ) L2(n) (7.4) 

:= hi +ii2- 
From the definition of G, 

h 1 = (H\v 1 -v 2 ) t , v'-v 2 )^ 

= (Jjci J xi (v 1 - v 2 )t, v 1 - v 2 ) L2(n) 

= {Jxi(v l - v 2 ) t , J x i(v l - v 2 )) L2{n) 



= {Jx^v 1 - v 2 ), Jxiiv 1 - v 2 )) L2(n) 



(7.5) 



L 2 (Q) 



Therefore, 



1 d 



h 1 >--\\j x ^v l -v 2 )\\ Z L2(n) -C(T,K) snp \\v 1 (s)-v 2 (s)\\ 2 L2{n) . (7.6) 

On the other hand, by Lemma 3.5, 

\Ii2\ < C\\g\j - gfjh^wWv 1 - v 2 \\ L 2 {n) \\vf\\ L 2 {n) 

< C(T, K) ( sup (\l\s) - l 2 (s)\ + \u\s) - W 2 ( S )|) ] sup \\v\s) - v 2 (s)\\ L 2 (n) 

\se[o,t] J se[o,t] 

< C sup Wv 1 ^) -v 2 {s)f L 2 m . 

se[o,t] 

(7.7) 
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I3 can be estimated similarly, 



\I 3 \ < Csup \\v\s) - v 2 (s)\\ 2 L 2 (R3) . 

[0,t] 



Therefore, 



^ll-W* 1 ~ v 2 )\\h ( n) + m^L 1 - L 2 \ 2 + ^[(JiR 1 - R 2 )) ■ (R 1 - R 2 )] 
< Csup\\v 1 (s)-v 2 (s)\\ 2 L 2 

[0,t] 

It follows that 

\\Jxiiv 1 - v 2 )(t)\\ 2 L2(n) + m\L x - L 2 \ 2 + {J{R l - R 2 )) ■ (R x - R 2 ) 

<C f sup ||^ 1 (s) - v 2 (s)\\ 2 L2m dT 
Jo se[0,r] 

Since X 1 is always a diffeomorphism and J is positive definite, thus 

ll^ 1 " Ahtw^) < C / sup ||^ - v 2 \\ 2 L 2 m (T)ds (7.8) 

By the Gronwall's inequality, v 1 = v 2 a.e. in [0, T] x R 3 . Uniqueness, as in [18] or 
earlier publication [21] for the Euler equations, combining the fact that the system is 
reversible, implies 

w eC([o,r];F s )nc 1 ([o,r];^ 1 ). 

In fact, from the preceding estimates, one can easily get that 

J t W v W 2 H°{n) + m J t \ L \ 2 + J t i jR - R ]^ c \Mh. {W Vv h~>(n) + IMIl2(r3) + l) . (7.9) 

It implies that once ||Vu||£,c»(m does not blow up, \\v\\h s will not blow up. Using the 
argument as in the paper [18] . one can get that the lifespan of the solution does not depend 
on s. 

Acknowledgements 

Gratitude is expressed specially to supervisor Professor Zhouping Xin and the referees 
for their careful reading of the manuscript and their fruitful suggestions. 



References 

[1] J. P. Bourguignon and H. Brezis, Remarks on the Euler equation, J. Fund. Anal. 
15(1974), 341-363. 



31 



[2] P. Cusmille and T. Takahashi, Wellposedness for the system modelling the motion of 
a rigid body of arbitrary form in an incompressible viscous fluid. Czechoslovak Math. 
J., 58(133)(2008), no.4, 961-992. 

[3] M. Dashti and J. C. Robinson, The motion of a fluid-rigid body system at the zero 
limit of the rigid body radius, Arch. Rational. Mech. Anal, 200(2011), no. 1, 285-312. 

[4] B. Desjardins and M. J. Esteban, Existence of weak solutions for the motion of rigid 
bodies in a viscous fluid, Arch. Ration. Mech. Anal, 146(1999), no. 1, 59-71. 

[5] B. Desjardins and M. J. Esteban, On weak solutions for fluid-rigid structure interac- 
tion: compressible and incompressible models, Comm. Partial Differential Equations, 
25(2000), no. 7-8, 1399-1413. 

[6] E. Dintelmann, M. Geissert and M. Hieber, Strong L p -solutions to the Navier-Stokes 
flow past moving obstacles: the case of several obstacles and time-dependent velocity, 
Trans. Amer. Math. Soc, 361(2009), no. 2, 653-669. 

[7] G. P. Galdi, An introduction to the mathematical theory of Navier-Stokes equations, 
Vol. I, Linearized steady problems, Springer and Verlag, 38(1994). 

[8] M. D. Gunzburger, Lee Hyung-Chun and G. Seregin, Global existence of weak solu- 
tions for viscous incompressible flows around a moving rigid body in three dimensions, 
J. Math. Fluid Mech., 2 (2000), no. 3, 219-266. 

[9] K. H. Hoffmann and V. N. Starovoitov, On a motion of a solid body in a viscous 
fluid. Two dimensional case, Adv. Math. Sci. Appl., 9(1999), no. 2, 633-648. 

[10] K. H. Hoffmann and V. N. Starovoitov, Zur Bewegung einer Kugel in einer zahen 
Flussigkeit[On the motion of a sphere in a viscous fluid], Doc. Math., 5(2000), 15-21. 

[11] J. G. Houot, J. S. Martin and M. Tucsnak, Existence of solutions for the equations 
modeling the motion of rigid bodies in an ideal fuild, J. Fund. Anal, 259(2010), 
2856-2885. 

[12] A. Inoue and M. Wakimoto, On existence of the solutions of the Navier-Stokes equa- 
tion in a time dependent domain, J. Fac. Sci. Univ. Tokyo Sect. IA Math., 24(1977), 
303-319. 

[13] N. V. Judakov, The solvability of the problem of the motion of a rigid body in a 
viscous incompressible fluid. (Russian) Dinamika Splosn. Sredy Vyp. 18 Dinamika 
Zidkost. so Svobod.Granicami (1974), 249-253, 255. 



32 



[14] T. Kato and C. Lai, Nonlinear evolution equations and the Euler flow, J. Fund. Anal. 
56(1984), no. 1, 15-28. 

[15] H. Kozono, On existence and uniqueness of a global classical solution of 2D Euler 
equations in a time dependent domain, J. Differential Eq., 57(1985), 275-302. 

[16] J. H. Ortega, L. Rosier and T. Takahashi, Classical solutions for the equations mod- 
elling the motion of a ball in a bidimensional incompressible perfect fluid, ESAIM: 
M2AN, 39(1) (2005), 79-108. 

[17] J. H. Ortega, L. Rosier and T. Takahashi, On the motion of a rigid body immersed in 
a bidimensional incompressible perfect fluid, Ann. I. H. Poincare, 24(2007), 139-165. 

[18] C. Rosier and L. Rosier, Smooth solutions for the motion of a ball in an incompressible 
perfect fluid, J. Funct. Anal, 256(2009), 1618-1641. 

[19] H. J. San Martin, V. N. Starovoitov and M. Tucsnak, Global weak solutions for the 
two dimensional motion of several rigid bodies in an incompressible viscous fluid, 
Arch. Rational Mech. Anal, 161(2002), no.2, 113-147. 

[20] D. Serre, Chute libre d'un solide dans un fluide visqueux incompressible. Existence, 
Japan J. Appl. Math., 4(1987), no. 1, 99-110. 

[21] J. Simon, Compact sets in the space L p (0,T;B), Ann. Mat. Pure Appl, 146(1987), 
65-96. 

[22] V. A. Solonnikov, Estimates in L p of solutions of elliptic and parabolic systems, 
Boundary Value Problems of Mathematical Physics. V. 1970, 157-185. 

[23] T. Takahashi and M. Tucsnak, Global strong solutions for the two dimensional motion 
of an infinite cylinder in a viscous fluid, J. Math. Fluid Mech., 6(2004), 53-77. 

[24] R. Temam, On the Euler equation of incompressible perfect fluids, J. Funct. Anal, 
20(1)(1975), 32-43. 

[25] Y. Wang and Z. P. Xin, On a weak solution for a two-dimensional fluid-rigid body 
system, submitted. 

[26] Y. Wang and Z. P. Xin, Analyticity of the semigroup associated with the fluid-rigid 
body problem and local existence of strong solutions, J. Funct. Anal, 261(2011), no. 
9, 2587-2616. 



33 



